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FOREWORD 


This  report  has  been  prepared  to  present  a  simple  modification  that  can  be  made  to 
the  Brouwer-Lyddane  artificial  satellite  theory  to  account  for  the  effects  of  atmospheric 
drag.  This  analytic  method  for  generating  ephemerides  for  satellites  has  been  implemented 
into  a  station  alert  prediction  algorithm  and  its  accuracy  determined  for  several  representative 
low-altitude  earth  satellites.  The  results  of  this  study  indicate  that  significant  improvement 
in  predicted  satellite  ephemerides  can  be  made  by  introducing  this  simple  drag  modification 
into  the  Brouwer-Lyddane  theory.  This  document  has  been  reviewed  and  approved  by  R.  J. 
Anderle. 


Released  by: 


J-.  /) /w 

O.  F.  BRAXTON,  Head 
Strategic  Systems  Department 
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INTRODUCTION 


The  Brouwer-Lyddane  theory  is  a  general  perturbation  method  that  can  be  used  to 
rapidly  predict  satellite  ephemerides.  This  theory,  however,  considers  only  terrestrial  oblateness 
perturbations  in  its  formulation  and  takes  no  account  of  atmospheric  drag  accelerations.  As  a 
result,  the  errors  induced  by  the  drag-free  Brouwer-Lyddane  theory  into  the  predicted  positions 
and  velocities  of  low-altitude  earth  satellites  can  be  quite  large  and  the  associated  ephemerides 
unacceptable  for  scheduling  ground  station-satellite  inview  periods.  The  inability  of  the  Brouwer- 
Lyddane  theory  to  provide  sufficiently  accurate  long-term  tracking  schedules  for  low-altitude 
satellites  can  result  in  undesirable  operational  difficulties. 

This  report  has  been  prepared  to  briefly  discuss  the  effectiveness  of  a  semi-empirical 
modification  that  can  be  made  to  the  Brouwer-Lyddane  theory  to  account  for  atmospheric 
drag  in  a  very  simple  way.  This  modified  theory  can  be  used  to  predict  station  alert  schedules 
for  low  altitude  earth  satellites  with  a  much  more  palatable  level  of  accuracy. 

A  discussion  of  the  formulation  of  the  semi-empirical  orbit  prediction  model  mentioned 
above  and  a  review  of  the  apparent  effectiveness  of  this  approach  in  diminishing  drag-induced 
errors  in  the  predicted  satellite  rise  time  that  initiate  a  station  inview  period  are  given  in  the 
following  sections. 


FORMULATION  OF  THE  DRAG-AUGMENTED  BROUWER-LYDDANE1  2  METHOD 


The  equations  given  below  are  those  of  the  drag-augmented  Brouwer-Lyddane  method 
and  consist  of  nothing  more  than  the  addition  of  time  rates  of  change  terms  to  the  original 
Brouwer  equations  for  the  osculating  semi-major  axis  (a)  and  eccentricity  (e)  and  for  the 
secular  mean  anomaly  (E").  The  remaining  equations  for  inclination  (i"),  argument  of  perigee 
(g)  and  right  ascension  of  the  ascending  node  (h)  are  unchanged.  These  additive  terms  in  the 
equations  below  are  enclosed  in  boxes  for  easy  identification. 


'Brouwer,  D.,  "Solution  of  the  Problem  of  Artificial  Satellite  Theory  Without  Drag,”  The  Astronomical  Journal,  Vol.  64,  No.  I  274. 
1 959,  pp.  378-397. 

2Lyddane,  R.  H.,  “Small  Eccentricities  or  Inclinations  in  the  Brouwer  Theory  of  Artificial  Satellites,"  The  Astronomical  Journal,  Vol. 
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Prior  to  the  transcription  of  the  drag-augmented  equations,  the  following  quantities 
must  be  defined: 


a" 

=  semi-major  axis  decay  rate 

e" 

=  eccentricity  decay  rate 

n 
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where  C/0 

(j  =  2,  3,  4,  5)  are  the  zonal 

harmonic  gravitational  expansion 

coefficients.  Then 

the  secular  terms  are  computed  from: 
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the  eccentric  anomaly  E'  is  obtained  from  a  Newton-Raplison  iteration  upon  the  Kepler  equation 
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The  final  osculating  values  for  a,  i.  and  li  are  computed  from  equations  (8).  (10).  and  (12), 
respectively.  Equations  (2).  (9),  (II),  and  (13)  are  used  to  calculate  final  osculating  values  for  8,  g, 
ind  e  from  the  following  relations: 
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e  sinfi"  +  e$  12  cos£"  , 
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and 
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All  of  the  input  data  required  to  generate  the  osculating  orbital  elements  from  the  above- 
augmented  equations;  i.e.,  a”,  a",  e"  e”,  i",  £",  li,  g^',  and  h",  can  be  obtained  directly  from  the 
NAVSPASUR  five-card  element  format  data.3  However,  the  semi-major  axis  and  associated  decay 
rate  provided  there  are  the  Kaula  semi-major  axis  (ak )  and  decay  rate  (ak )  expressed  in  earth  radii. 
The  Kaula  semi-major  axis  can  be  converted  to  the  desired  linear  units,  when  one  uses  the  relation: 
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(22) 


where  ae  is  the  earth’s  semi-major  axis, 
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and  J2  =  •  zonal  harmonic  gravitational  constant. 

Similarly,  ti.  esired  decay  rate  a"  can  be  obtained  from  aR  by  taking  the  time  derivative  of 
equation  (22).  This  ves 
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^McCandless,  A.  R.,  “NAVSPASUR.  SPADATS  Type,  5-Card  Element  Format;  description  of,”  memorandum  V22 1  :F.AT:jmd, 
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where 


X  = 


(25) 


In  the  last  expression  it  has  been  assumed  that 
(1")  =  0. 


(26) 


EVALUATION  METHOD 


In  order  to  evaluate  the  utility  of  the  drag-augmentation  method,  five-card  element 
sets  were  obtained  from  NAVSPASUR  for  two  92-min.  period  satellites  at  7-day  intervals 
for  the  timespan  from  9  January  1981  through  24  April  1982.  During  this  timespan,  these 
two  satellites  had  perigee  altitudes  which  remained  in  a  160-nmi  to  180-nmi  altitude  box. 

These  data  were  used  to  generate  30-day  predicted  station  alert  schedules  every  7  days 
during  the  timespan  for  a  single  tracking  station.  The  satellite  rise  times  generated  near  the 
epochs  of  the  7-day  element  sets  were  assumed  to  be  the  true  rise  times,  tRx.  These  were 
compared  with  the  associated  predicted  rise  times,  tR|>.  generated  for  30  days  at  7-day  epoch 
intervals.  Orbit  adjust  schedules,  obtained  for  these  two  satellites  for  the  span  of  interest, 
were  used  to  eliminate  those  values  of  tRp  which  were  invalidated  due  to  orbit  adjusts 
occurring  at  some  time  during  the  prediction  interval.  The  remaining  valid  tRp  were  differenced 
with  the  associated  tR  T  and  their  absolute  values  averaged  as  a  function  of  prediction  interval. 


RESULTS  AND  CONCLUSIONS 


The  results  of  the  averaging  procedure  described  in  the  previous  section  are  shown  in 
Figure  1  along  with  the  associated  2o  confidence  envelope.  Also  included  on  this  figure  is  a 
plot  of  the  errors  obtained  from  a  typical  drag-free  Brouwer-Lyddane  prediction  for  one  of 
the  two  satellites  discussed  previously  when  its  perigee  altitude  was  ~  1 78  nmi  and  its  argument 
of  perigee  was  on  the  sunlit  side  of  the  earth.  It  is  apparent  from  these  three  curves  that 
the  inclusion  of  the  simple  drag  decay  terms  in  the  Brouwer-Lyddane  theory  has  dramatically 
improved  the  long-term  quality  of  the  predicted  satellite  rise  times. 
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